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I.  Introduction  and  Summary. 


The  present  report  is  a  continuation  of  an  earlier  work  involving 
some  Issues  on  a  Cs-Ne  high  power  flashlamp  and  its  modeling.1  Because  of  the 
complexity  of  the  physics  encountered  in  an  arc  discharge  of  high  gas  pressures 
and  low  voltages,  it  is  necessary  to  divide  the  theoretical  effort  into  two 
parts:  microscopic  and  macroscopic  processes.  A  detailed  discussion  of  this 
modeling  has  been  given  separately  in  a  status  report. 

>!/ 

In  the  present  report,  we  are  primarily  concerned  with  atomic 
processes  in  a  high  pressure  gas  discharge.  The  theoretical  endeavour  is 
highlighted  in  a  detailed  calculation  of  spectral  line  profiles,  which  are 
essential  to  meaningfully  model  an  arc  lamp  performance.  We  also  examine 
the  effect  of  neon  as  a  buffer  gas  in  a  high  pressure  cesium  discharge  tube. 
Under  the  condition  of  comparable  cesium  and  neon  pressures,  we  concluded 
that  neon  buffer  gas  Is  beneficial  to  the  overall  lamp  performance.  The  de¬ 
tailed  discussion  can  be  found  In  Reference  2. 
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This  report  is  devoted  to  a  rigorous  analysis  and  calculations 
of  spectral  line  shapes  for  high  pressure  and  low  voltage  arc  lamp  applications. 
In  Section  II,  a  theoretical  description  of  spectral  line  shape  analysis  is 
given.  The  formultion  used  is  a  straightforward  generalization  of  those  used  in 
a  previous  report.1  Section  III  describes  spectral  line  shape  analysis  and 
calculations.  The  calculations  include  the  following  four  types  of  interaction: 

•  neutral -neutral  interactions:  Cs-Cs  and  Cs-Ne, 

•  neutral -charged  Interactions:  Cs-e  and  Cs-Cs+. 

These  collisions  lead  to  broadenings,  shifts,  and  asymmetries  of  the  observed 
spectral  line  shapes  and  to  the  occasional  appearance  of  diffuse  bands. 

In  Section  IV,  our  future  research  plans  and  recommendations 
are  indicated.  Finally,  a  listing  of  the  code  developed  for  applying  the  given 
model  is  presented  in  Appendix  A  followed  by  a  detailed  mathematical  derivation 


of  spectral  line  shapes  in  Appendix  B*. 

Based  on  our  limited  efforts,  we  summarize  our  theoretical 
findings  as  follows: 

•  for  high  power  Cs-Ne  flash! amps,  the  neon  pressure 
should  be  comparable  to  cesium  pressure, 

•  if  neon  (or  other  inert  gases)  pressure  is  several  orders 
of  magnitude  higher  than  cesium  (or  other  alkali  atoms) 
pressures,  an  infrared  laser  seems  possible  by  optical 
pumping  methods,  and 

•  for  high  power  Cs-Ne  arc  lamps,  the  spectral  profile  is 
mainly  determined  by  cesium  doublet  mixing  effect. 


II.  Theoretical  Description. 


The  study  of  spectral  line  shapes  was  actively  pursued  during 
the  past  fifty  years.  Numerous  review  articles  and  books  have  been  written 
since  late  1930.  **  Regardless  of  theoretical  approaches  taken  in  the  pre¬ 
vious  work,  it  is  assumed  that  an  Interaction  potential  between  a  radiating 
particle  (atom  or  molecule)  and  a  perturbing  particle  as  a  function  of 
interparticle  distance  is  known.  As  a  result,  experimental  measurements  of 
spectral  line  shapes  are  used  to  determine  the  parameters  introduced  in  the 
assumed  interaction  potential.  This  method,  however,  achieved  only  a  limited 
success.  Several  attempts  were  made  in  recent  years  to  fit  experimental  line 
shapes  by  introducing  more  and  more  adjustable  parameters  into  the  existing 
thoery. 

The  effort  of  this  kind  is,  however,  fruitless  in  our  opinion. 
The  reason  Is  simple.  First  of  all,  the  interaction  potential  cannot,  in 
general,  be  represented  as  a  function  of  particle  separations  in  a  simple 
mathematical  form.^  Secondly,  from  a  rigorous  theoretical  point  of  view 


* 


Appendix  B  not  reproduced  in  all  copies  of  this  report. 


It  would  be  highly  unsatisfactory  because  of  its  lack  of  predictability  if 
it  should  contain  adjustable  parameters. 


In  order  to  circumvent  the  deficiency  as  stated  above,  a  new 
approach  must  be  developed.  The  method  of  calculating  spectral  line  shapes 
to  be  given  below  has  two  important  advantages  over  the  previous  approach. 
Firstly,  no  assumed  form  of  interaction  potential  as  a  function  of  inter¬ 
particle  separation  is  used.  Secondly,  the  formulation  can  be  easily  extended 
to  interactions  involving  various  different  species  as  Ions  as  the  interaction 
is  of  electromagnetic  origin.  As  a  consequence,  the  resulting  line  shape  does 
not  contain  any  adjustable  parameters.  It  is  completely  determined  by  the  atomic 
or  molecular  structures  of  the  radiating  and  perturbing  particles. 

Calculations  of  spectral  line  shapes  can  be  accomplished  by  a 
simple  concept.  For  an  atom  or  a  molecule  initially  in  a  given  excited  state 
making  a  radiative  transition  to  a  lower  quantum  state,  a  photon  is  emitted. 

The  frequency  distribution  of  the  emitted  photon  for  a  given  transition  is 
called  a  spectral  line  shape.  For  a  single  Isolated  atom  without  any  external 
perturbation,  the  line  shape  can  be  derived  using  the  well-known  method  of 
Weisskopf  and  Wigner.6  This  method  can  be  generalized  to  the  problem  including 
the  external  perturbations  due  to  interaction  with  other  atoms  of  the  same  or 
different  species.  The  effect  of  external  perturbation  is  to  modify  the  upper 
and  lower  radiating  states.  In  the  following  section,  we  indicate  how  to  carry 
out  an  explicit  calculation  of  spectral  line  shapes. 
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III.  Spectral  Line  Shape  Analysis  and  Calculations 


We  shall  consider  the  far-infrared  radiation  from  the  cesium 
atom  in  the  wavelength  range  between  lp  and  5p.  These  radiations  involve 
7S  -*•  6P,  7P  -*•  5D  and  5D  -*•  6P  as  transitions.  These  transitions 
do  not  involve  the  ground  state  6S  as  the  lower  level.  Thus  the  photon  dis¬ 
persion  is  likely  to  be  unimportant.  Assuming  this  to  be  the  case,  we  may 
approximate 


4*.  *  |k(  =  w. 

k 


As  a  consequence,  the  photon  wavefunction  renormalization  constant  Zg(u>  )  is 
equal  to  unity,  and  Equation  (6 f  becomes,  using  the  nonrelativistic  approximation. 


*  e  (2“V>'U2  <n£  •  *A  • 


\III.1) 


The  electron-atom  interaction  kernel  KaD  in  the  same  approximation 


can  be  written,  neglecting  P^(q ,q0) 


Kn, n«  (p’q)  "  *  V  Ve  “l|n2>/iq  -Q0  )  .  * 

1  2  qo=i  h  ,n2  +  7  V 


(III- 2) 


The  atom-atom  interaction  kernel  Km  can  likewise  be  obtained. 


K^,n2,„3,n4'2-5)  *  -  f  <»il  <n3i\e 


n,,n?  2  -1 

x  [ 2)  -  Q2] 


(III. 3) 


If  now  we  use  the  dipole  approximations.  Equations  (III. 1-1 I I .3)  take  the  form 


*  Equations  ( 1 ) - ( 90 )  refer  to  Equations  in  Appendix  B. 
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where  d„  _  is  the  dipole  moment  defined  as 

nj,n2 


d  - -  <n,  |r|n9>  . 

nrn2  ^  12 


We  shall  use  the  above  approximation  to  analyze  the  quantities  appearing 
in  our  lineshape  formula.  The  quantity  Ye  In  Equation  (64)  is  due  to  the  photon- 
atom  interaction,  which  is  obtained  from  Equation  (33).  Explicitly, 


=  2  I  ;ao  I'3!*.  !2  e  (E?  -  O 
e  3  1  e,n 1  e,n 1  e  n' 


which  is  a  sum  of  Einstein  A-coefficients  for  dipole  allowed  transitions  to  the 
lower  states. 

The  quantities  that  appear  in  Equations  (69)  and  (73)  are  due  to  electron- 
atom  and  atom-atom  Interactions.  From  Equations  (73),  letting  n  -  0,  we  have 
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TJ,(o) 
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(III. 4) 
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In  Equation  (III. 4),  g(p)  is  the  electron  distribution  function,  which  is 
assumed  to  be  Maxwellian  with  an  electron  temperature  T„  , 

P 

9<P>  »  (.  p |r3/Vf2/p«  ) 


PB  =  (2  m  kBT6) 


The  effect  of  electron  drift  velocity  will  be  ignored.  Both  functions  defined  in 
Equltion  (111,5)  are  complex  with  positive  imaginary  parts,  with  y^  being  the  radiative 
decay  width  of  the  radiating  state  |j>.  The  interaction  kernels  iK^fand  |Ka6|2 
are  functions  of  Q  and  q  respectively.  They  are  independent  of  P  and  P  .  Further- 
more,  f(p)  and  g(p)  are  also  independent  on  the  direction  of  P  and  p.  As  a  con¬ 
sequence,  the  angular  integrations  in  Equations  (III. 4)  and  (III. 5)  can  be  carried  out. 
Let 


A(P  )  *  2u  /d3P  f (P)f  -r 

0  J  J  (  2tt)  3  0Z+ 


(2tt)j  Q^+2(^-f^0)  •  ff  +a 


where 


Im  a  <  0 


Carrying  out  the  angular  integrations,  we  get 
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If  f(P)  is  assumed  to  be  Maxwellian  as  in  Equation  (79),  then  Equation  (III. 6) 
the  form 


A(Pq) 


-P  /P 


(2tt)5/2P  p  J 

o  a 


a  P+P. 


TP+PTT  B(lp+p0|) 


(III. 7) 


where 


00 

B(x)  =  /  QdQ | K(Q) | 2  In  (Q2+2Qx  +  a) 


(III. 8) 


Consider  the  electron-atom  interaction.  If  the  dipole  approximation  is 
made,  then 


|K(Q)!2  oc  Q~2 


The  integrand  in  Equation  (III. 8)  behaves  as  f dr  large  Q, 


In  Q/Q 


and  the  integral  diverges.  The  divergence  of  this  kind  is  encountered  in  the  usual 
thoery  based  on  a  multipole  expansion  of  the  static  coulomb  potential,  and  adjustable 
parameters  are  often  invoked  to  make  the  results  finite. 


...  {Hu’jwf-'jw.'j  .umw «  ^  «j.  ^Aiv.vgvj,v;*  .va  *•  :.’•  -1 
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“  "  v i* 


In  the  present  formulation,  this  divergence  is  seen  from  the  approximation 
of  the  matrix  elements  of  the  form 


<n1|eiQz|n2>  s  i  Q  <ni|2  |n2> 

This  approximation  is  correct  only  when  Q  is  small.  For  large  Q,  the  matrix 
element  must  decrease  at  least  as  Q"1.accordino  to  the  Riemann-Lebesgue  lermw.  Because 
of  this  fact,  we  showed  in  an  earlier  report  that  we  may  approximate 

2 

<n1|ei^z|n2>  s  i  Q  <nj|<|n2>/[T  +  f-Untlz|n2>|2]  (III. 9) 

If  Equation  (III. 9)  is  used  to  evaluate  electron  excitation  cross  sections,  an 
agreement  with  experimental  data  is  achieved.  Note  that  Equation  (II  1.9)  is 
consistent  with  the  expected  Independence  for  both  large  and  small  Q. 

If  now  we  used  Equation  (III. 9)  to  evaluate  |K(Q))|2,  than 


|K(Q)|2-  ,  1  ,,  (HI. 10) 

<r(i  +  birr 

The  integral  in  Equation  (III. 8)  thus  converges  rapidly  for  large  Q.  Note  that 

there  is  a  logarithmic  singularity  in  Equation  (III. 8)  at  Q  *  0  if  Equation  (III. 10) 

is  used.  This  arises  from  the  Q’2  dependence  in  the  electron-atom  interaction  kernel 

where  the  lifetime  of  the  coupled  radiating  state  |j>  has  been  ignored.  If  this 

-2 

finite  lifetime  and  retardation  effect  are  taken  into  account,  then  Q  should  be 
replaced  by 

«*  -  <4j.n'  *  1  V2)'1 

and  the  logarithmic  singularity  on  the  real  axis  disappear*. 
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For  atom-atom  interaction,  we  have,  if  dipole  approximations  are  used. 


|K(Q) : 2  =  0  (1) 


for  large  Q.  Then  the  integral  in  Equation  (III. 8)  diverges  quadratically. 

However,  using  Equation  (III. 9)  and  taking  into  account  the  finite  life-time 
of  the  radiating  state,  a  well-defined  finite  physical  result  is  again  obtained. 

According  to  the  foregoing  discussion,  the  functions  K(Q)'2  in  Equation  (III. 8) 
have  the  following  forms: 

t  Electron-atom  interaction 


|K(q)|2  =  e4|<j|z|n'>|2q^  [(<HVjjn.)2  +  ^][(q+^  n«)2  +  -^] 


*  (1  +  ffNUIn'  , V)*1 


(III.  11) 


•  Atom-atom  interaction 


I K(Q) 1 2  *  e4|<j|z|n'><g|z|n">|2  G(Q) 

«((1*  C2„,Q2)2(1  ♦  02„Q2)2C(Q  '  „')2  *  J1  3 

Y  2 


(III. 12) 
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1.12),  averages  over  dipole  orientations  have  been  performed. 

present  treatment,  we  also  include  the  retardation  effect  as  in 
11)  and  (I II. 12).  This  effect  has  not  been  included  in  existing  theories 
eshape  as  far  as  we  know.  However,  as  pointed  out  by  Stevens,  the 
ect  is  important  in  determining  the  interaction  energies  between  two 
separation.  This  effect  should  lead  to  a  certain  modification  on 
arameters  in  particular,  and  the  lineshape  itself  in  general. 

e  calculation  of  spectral  line  profiles  can  be  performed  if  the 
ation  (67),  are  satisfied.  In  this  connection,  the  lineshape 
Dn  (68),  may  be  used.  From  Equations  (74)  and  (68),  may  be  used. 

(74)  and  (68),  after  carrying  out  the  angular  integration,  we 
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(III. 13) 
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jgrating  Equation  (II I. 13)  by  parts,  we  have 
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(III. 15) 


N±(P0.Z)  *  R2(P0)  j?-  R!<P0)  -  [z  *  R2(P0)  *  ar  Rz(po> 


PQ  dependence  in  R1(PQ)  and  R2(PQ)  is  ignored.  Equation  (III. 14)  reduces  to  the 
1- known  Voigt  profile. 


To  evaluate  the  integrals  in  Equation  < I I I . 14)  as  a  function  of  Z,  it  is 
;ssary  to  calculate  Rj^q)  and  R2(PQ^’  W^1C^  are  expressed  in  terms  of  Ag 
An  as  defined  in  Equation  (69).  They  are  given  explicitly  below. 

Define  the  reduced  masses  jj  and  u„  for  atom  and  electron  respectively, 
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We  write 
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The  functions  B^(|P  +  PQ|)  and  B^(|p  +  PQ|)  are  computed  according  to  Equation  (III. 8) 
with  the  parameters  occuring  there  identified  as  follows:  For  Bg,  use  Equation  (III. 11) 
for  |K(Q) |2,  and 


ls¥n',j*!V‘VilP( 


For  b£,  use  Equation  (III.  12)  for  |K(Q)|2,  and 


a 


i 

1 


The  computer  program  is  given  in  Appendix  A.  Certain  sample  calculations 
are  shown  in  Reference  2. 
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IV.  Future  Plans 


For  a  single  isolated  line,  a  general  expression  of  spectral  line 
shapes  was  obtained,  as  shown  in  Equation  (72)  of  Appendix  B?  Because  of  the 
limited  funds  available  at  present,  a  general  numerical  evaluation  of  Equation  (72) 
cannot  be  performed.  Only  the  approxi rated  expression.  Equation  (68),  to  the 
exact  formula.  Equation  (72),  was  numerically  implemented.  There  are  circum¬ 
stances  in  which  the  approximations,  Equations  (67),  leading  to  Equation  (68) 
may  fail..  In  this  connection,  numerical  calculations  of  the  exact  expression. 
Equation  (72),  must  be  carried  out. 

In  addition  to  the  above  technical  program  aimed  at  performing  an 
accurate  spectral  line  shape  calculation,  there  are  three  issues  needed  to  be 
addressed  for  future  program  needs.  They  are 

A.  Physics  Issues 

•  Atomic  Model  Development 

•  Atomic  Core  Polarization  Effects 

•  Spectral  Line  Profile  Formation 

•  Electrical  Conducti  vity  Calculation 

•  Thermal  Conductivity  Calculation 

B.  Devise  Improvement  Issues 

•  Emission  Spectrum  Tailoring 

•  Flashlamp  Reliability  Improvement 

•  Optical  Gas  Mixture 

•  Cesium-Neon  Discharge  Laser 

•  Cesium  Vapor  Heating  Problems 


*  Not  reproduced  in  all  copies  of  this  report. 
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Diagnostic  Issues 


•  Discharge  Plasma  Diagnosis 

•  EO  Systems  Quality  Control 


Detailed  discussions  of  the  above  issues  can  be  found  in  Reference  2. 


Appendix  A 


Computer  Program  Description  and  Listing 

A. 1  Program  Description 

This  is  a  description  of  the  computer  routines  mentioned  in 
in  Section  III  above  which  are  used  to  compute  line  shape  by 
solving  equation  III. 13.  The  description  of  the  calculations  in 
each  routine  follows  the  order  of  the  calculations  in  that  routine. 
The  routines  appear  in  the  order  that  they  are  called  in  the 
program.  A  flowchart  for  subroutine  SHPE  and  its  associated  routines 
follows  in  Figure  1* 

Subroutine  SHPE  is  the  primary  routine  -  it  calculates  the  final 
lineshape: 

Subroutines  INTLZ.SETDAT  and  MED  are  called  to  set  the  physical 
constants  and  atomic  data. 

For  each  allowed  transition  between  fine  structure  lines  a 
lineshape  is  calculated: 

Subroutine  SETTRN  is  called  to  set  indices  corresponding  to 
to  the  fine  structure  lineshape  being  calculated. 

Subroutine  RS  is  called  to  provide  the  values  of  the 
parameters  Rl(Po)  ,R2(Po) ,  (from  equation  III. 13),  at  the 
Hermite  abscissa  points. 

At  each  freouency  point  used  to  define  the  lineshape: 


The  integral  in  equation  111.13  is  performed  by  Hermite 
integration  to  produce  a  lineshape  value. 


The  weighted  value  is  added  to  the  final  lineshape  at 
that  frequency  point. 

Subroutine  RS  calculates  the  values  of  the  parameters  Rl(Po),R2(Po) 
at  the  Hermite  abscissa  points: 

First  the  contribution  to  Rl(Po),  R2(Po)  from  the  initial  atomic 
state: 

Subroutine  SET  is  called  to  set  parameter  values  corresp- 
pondina  to  which  atomic  state  and  exciting  states  are 
invol ved. 

Subroutine  AS  is  called  for  each  of  the  species  (atoms, 
electrons,  buffer  atoms,  ions)  to  provide  the  values 
of  the  parameters  Ae,An  as  defined  in  equation  (69). 

Sun  the  Ae,An  parameter  contribution  from  each  species. 

Then  the  contribution  to  Rl(Po) ,R2(Po)  from  the  final 
atomic  state:  repeating  the  above  procedure. 

The  values  of  the  parameters  Rl(Po),R2(Po)  for  the  initial 
and  final  states  are  summed. 

Subroutine  AS  calculates  the  value  of  A^  as  shown  in  the  integrals 
on  page  12  by  performing  a  Hermite  integration. 

(*  -  defines  what  species  is  involved.) 


Subroutine  SUMRI  is  called  separately  to  evaluate  the 
imaginary  and  real  parts  of  the  function  B  at  Hermite 
abscissa  points. 

Subroutine  Herm  is  called  to  perform  the  integration. 

Subroutine  SUMRI  calculates  the  real  or  imaginary  part  of  the 
function  B  by  evaluating  the  integrals  as  shown  in  equation  III. 8  . 
It  performs  a  trapaziodal  method  of  integration: 

First  it  breaks  the  range  of  integration  into  one  hundred  steps. 

At  each  point  subroutine  NBEQ  is  called  to  evaluate  the  function 
that  is  being  integrated  and  the  Integral  to  that 
point  is  summed. 

Subroutine  NBEQ  evaluates  the  function  that  is  being  integrated  in 
equation  III. 8  at  a  specified  integration  point. 

First  it  determines  if  its  the  real  or  imaginary  part  of  the 
calculation  since  different  formulas  are  required  and  then 
procedes  to  the  apropriate  part  of  the  subroutine. 

Subroutine  KS  is  called  to  evaluate  the  appropriate  value  or 
I  K(Q)  j2  as  shown  in  equations  III. 11,  III. 12. 

p 

The  appropriate  value  of  J K( Q ) |  is  used  to  calculate  B  . 

Subroutine  KS  calculates  the  value  of |K( Q) | ^  : 

2 

It  determines  what  formulation  for  |K(Q)|  should  be  used. 

For  a  charged  particle  equation  III. 11  Is  solved. 

For  a  neutral  particle  equation  III. 12  is  solved. 


INTLZ 
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A. 2  Computer  Program  listing 


SUBROUTINE  SHPE(ZRNGE» SHAPE 'RHAtRHB'RHGfRHD? 

1  TEMA»TEMB»TEMG»TEMD»FO»IS»NZ) 

IMPLICIT  REAL*8(A-HrO-Z) 

COMMON  /RSVAL/  Rl(20)fR2(20> 

COMMON  /HERDAT/  HWT(20>»HABS(20)»NHERM 
COMMON  /TRANS/  INIT»INITXiLAST»LASTX»DJE»DJNt 
1  6IN»GAME»SHIFTZ»NTRANS(30) rNEXCIT(5) » FRAC? IDK 

COMMON  /CSNE/  AMASS r EKASS 1 6HASS iZETAr RL J 
COMMON  /NEDM/  RH0A»RH0BiRH06rRH0D» 

1  TA»TB»TGrTD» 

2  PA»PB?PG?PD» 

3  RKA»RKB»RKG»RKD 
COMMON  /CONST/  E2» PI fVELC*RVKB» CONST 1 
DIMENSION  ZRNGE (NZ)» SHAPE (NZ) 

DIMENSION  FI (100) 

SUBROUTINE  SHPE  -  DRIVING  ROUTINE  TO  CALCULATE  LINESHAPE 
CALL  INPUTS 

ZRNGE (NZ)  -  ARRAY  OF  FREQUENCIES  (  ) 

RH*  -  SPECIES  DENSITY  (  /CM3) 

TEN*  -  SPECIES  TEMPERATURE  (DEG  K> 

(t  -  A(ATON)»  B( ELECTRON)?  G( BUFFER)?  D(ION)  ) 

FO  -  BOHR  FREQUENCY  (  /CM) 

IS  -  INDEX  TO  IDENTIFY  TRANSITION 
NZ  -  NUMBER  OF  FREQUENCIES 

CALL  OUTPUTS 

SHAPE (NZ)  -  LINESHAPE 

VARIABLES 

ZETA  -  DOPPLER  WIDTH  (  /CM) 

NTRANS(IS)  -  NUMBER  OF  FINE  STRUCTURE  LINES 

DATA  KALL/O/ 

IF(  KALL  »EQ.  1  )  60  TO  20 
KALLS1 
CALL  INTLZ 
CALL  SETDAT 
20  CONTINUE 
RHOA*RHA 


i 
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RHOG-RHG 

RHOD*RHD 

TA*TEMA 

TB*TEMB 

TG«TEMG 

TD*TEMD 

CALL  NED 

ZET A*P A*2 .0*PI*F0/< VELCtAKASS ) 

NTRAN=NTRANS(IS) 

DO  50  IZ«1>NZ 
SHAPEUZ>«0.0 
SO  CONTINUE 
C 

C  LOOP  OVER  LINE  TRANSITIONS 
C 

DO  100  NT*1»NTRAN 
CALL  SETTRN(NTrlS) 

CALL  RS 

PRINT  1010* <R1 (II) »R2(II)» 11*11 »20) 

1010  FORMAT ( 1X> 1P10E12 «2) 

C 

C  LOOP  OVER  FREQUENCY  VALUES 

C 

DO  200  IZ*1'NZ 
Z=ZRNGE ( IZ) -SHIFTZ 
F1<IZ)«0.0 
C 

C  HERNITE  INTEGRATION 

C 

DO  300  IH»ltNHERM 
W»HABS(IH) 

WT«HWT(IH) 

IF(M.LT.O.O)  GO  TO  300 

Fl(IZ>*Fl<IZ)+Mt<DATAN((Z+Rl(IH>+ZETA*H>/R2(IH)> 
1  -DAT AN ( (ZtRl(IH)-ZETAtU)/R2(IH) ) )*WT 

300  CONTINUE 

SHAPE(IZ)*SHAPE(IZ)+Fl(IZ)tFRAC 
200  CONTINUE 
100  CONTINUE 
RETURN 
END 

SUBROUTINE  INTLZ 
IMPLICIT  REALM (A-HiO-Z) 

COMMON  /CSNE/  AMASS  >  EMASS 1 6MASS  » ZET A » RL J 


uuuuuuuu  uuuuuuuuu 


COMMON  /CONST/  E2.PI tVELCfRVKBtCONSTl 
COMMON  /HEW/  RHOA  t  RHOB  t  RHOG  t  RHOD  t 

1  TAtTBtT6tT0t 

2  PAtPBtPGtPDr 

3  RKAtRKBtRKGtRKD 

COMMON  /TRANS/  INITtlNITXfLASTtLASTXtDJEfDJNt 
1  GINtGAMEtSHIFTZtNTRANS(50)tNEXCIT(5)tFRACtIDK 

ftt  EMASS  -  ELECTRON  MASS  (CM-1) 
tff  AMASS  -  ATOMIC  MASS  (CM-1) 
fft  GMASS  -  BUFFER  MASS  (CM-1) 
m  CONSTI- 

m  VELC  -  VELOCITY  LIGHT  (CM/S) 
tff  RVKB  -  BOLTZMANN  CONSTANT 


E2*l. 0/137.0 
EMASS-1. 0/3. 861 60-11 
AMASS=1 33 . Of 1836 . Of EMASS 
GMASS=20< f 1836. Of EMASS 
C0NST1*6.73E~14 
VELC* 3. OE 10 

RVKB=1 .38070-16/1.98650-16 

PI*3. 1415926500 

RETURN 

END 

SUBROUTINE  SETDAT 
IMPLICIT  REALt8(A-H»0-Z) 

COMMON  /TRANS/  INIT » INITXt LAST rLASTXtDJEr 0 JNi 
1  6IN»GAME»SHIFTZ»NTRANS(50) »NEXCIT(5) »FRAC» IOK 

COMMON  /ATOAT/  NSHELL(60).LM0(60)iJM0(60)»EN(6O).GAMARR(60) 
COMMON  /CONST/  E2>PIrVELCrRVKBrCQNSTl 

SUBROUTINE  SETDAT  -  SETS  ATOMIC  DATA 


NSHELL(I)  -  PRINCIPAL  OUANTUM  NUMBER 
LMO(I)  -  ORBITAL  ANGULAR  MOMENTUM 
JMO(I)  -  TOTAL  ANGULAR  MOMENTUM  f  2 
EN(I)  -  ENERGY  LEVEL  OF  STATE  (  /CM) 
GAMARR(I)  -  LIFETIME  OF  STATE  (  /S> 


DATA  NSHELL 

/  6r 

6» 

6* 

5t 

5t 

7t 

7 1 

7t 

1 

6» 

6» 

8  r 

4t 

At 

8t 

8t 

7 1 

1 

7t 

9t 

5» 

5? 

5t 

9t 

9t 

8t 

1 

8 » 

10. 

6» 

6t 

6t 

6t 

lOt 

10t 

1 

9» 

9i 

lit 

7f 

7t 

lit 

lit 

10t 
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DATA  LHO 


DATA  JKO 


DATA  EN 


DATA  6AHARR 


10? 

9? 

12?  8? 

9? 10*0/ 

8? 

12? 

12? 

11? 

11? 

0? 

1? 

1? 

2? 

2? 

0? 

1? 

1? 

2? 

2? 

0? 

3? 

3? 

1? 

1? 

2? 

2? 

0? 

3? 

3? 

4? 

1? 

1? 

2? 

2? 

0? 

3? 

3? 

4? 

5? 

1? 

1? 

2* 

2? 

0? 

3? 

3? 

1? 

1? 

2? 

2? 

3? 

0?  3? 

3? 10*0/ 

3? 

1? 

1? 

2? 

2? 

1? 

1? 

3? 

3? 

5? 

1? 

1? 

3? 

3? 

5? 

1? 

7? 

5? 

1? 

3? 

3? 

5? 

1? 

7? 

5? 

79? 

1? 

3? 

3? 

5? 

1? 

7? 

5? 

79? 

911? 

1? 

3? 

3? 

5? 

1? 

7? 

5? 

1? 

3? 

3? 

5? 

7? 

1?  7? 

5? 10*0/ 

5? 

1? 

3? 

3? 

5? 

/  0.00D0?11178.24D0? 11732. 3500? 14499. 4900? 

14597. 08D0?18535.51D0?21765.65D0?21946.66D0? 
22588 . 89D0 ? 22631 . 83D0  ?  2431 7 . 1 7DO ? 24472 . 29D0 ? 
24472 • 46D0 ? 25709 • 1 4D0 ? 25791 • 78D0 ? 26047 • 8600 ? 
26068 • 83D0 » 2691 0 . 6800 ? 24971 . 42D0 » 26971 . 56D0 ? 
27010. OODO? 27637 • 29D0>  27681 . 96D0 ? 2781 1 • 25D0  ? 
27822 . 94D0 » 28300 . 28D0 ? 28329 . 66D0 ? 28329 • 76D0  ? 
28347 . OODO ? 28356 • OODO » 28727 . 09D0 ? 28753 • 93D0  ? 
28828 • 90D0 ? 28836 • 06D0  > 291 30 . 00 DC » 291 48 • 1 6D0 ? 
29148. 23D0 » 29403 . 68D0 ? 29421 « 10D0 » 29468 . 54D0 » 
29473 . 22D0f 29666 . OODO  r 29678 . 94D0 > 29678 . 98D0 1 
29852 . 85D0f 29864 . 72D0  > 29896. 64D0» 29899 . 89D0 . 
30042 . 52D0 1 30042 . 54D0 » 10*0 . DO/ 

/  0.  t  3.280D+07.  3.7400407?  8.620D405? 

8.620D+05?  1.897D+07?  7.325D+06?  7.325D+06? 
1.7C0D+07?  1.700D+07?  1.044D+07?  1.630D+07? 
1.630D407?  2.601D+06?  2.601D406?  1.092D+07? 
1.092D407*  4 .6400406?  4.317D406?  4.317D406? 
2.630D406?  2.250D405?  2.250D405?  6.550D405? 
6.550D406?  2.457D406?  1.665D406?  1.665D406? 
4.870D405?  4.870D405?  3.731D406?  3.731D406? 
3.731D406?  3.731D406?  3.731D406?  7.434D405? 
7.434D405?  7.434D405?  7.434D405?  2.440D406? 
2.440D406?  2.440D406?  3.738D405?  3.738D405? 
3.738D405?  3.738D405?  8.130D405?  8.130D405? 

1 • 780D405  ?  1 .7800405? 1040./ 


DIVIDE  BY  VELOCITY  OF  LIGHT  -  GAMARR(I)  <  /CH) 


'  «**  »*■*  »  ' 


?•'  V* 


DO  200  N*1»A0 
GAMARR(N)*GAMARR<N)/VELC 
200  CONTINUE 
RETURN 
END 

SUBROUTINE  NED 
IMPLICIT  REAL*8(A-H»0-Z) 

COMMON  /CSNE/  AMASSrEMASSfGMASStZETAiRLJ 
COMMON  /ATDAT/  NSHELL(60)»LM0(60)»JN0(60)»EN(60)*GANARR(60) 
COMMON  /TRANS/  INITrINITX»LASTtLASTXfDJE*DJN» 

1  GIN»GAME»SHIFTZ»NTRANS(50) .NEXCIT(3) tFRACt IDK 

COMMON  /HEDM/  RHOA  t RHOB  t RHOG  t RKOD  > 

1  TAfTBtTGfTDt 

2  PA»PB*PG»PD» 

3  RKAiRKB'RKG'RKD 
COMMON  /BCON/  RITAA»RITAB»RITAG»RITAD» 

1  EPSA'EPSB'EPSG»EPSDf 

2  RKSIA»RKSIB»RKSIG»RKSID» 

3  CJA»CJB>CJGfCJD> 

4  6AHJfGt6P»DNPJ»DJ«LJ 
COMMON  /CONST/  E2 i P I ! VELC i R  VKB i CONST 1 

C  SUBROUTINE  MED  CALCULATES  PARAMETERS  UHICH  REMAIN  CONSTANT 
THROUGH  EACH  CALL  TO  SUBROUTINE  SHPE 

P *  -  SPECIES  MOMENTUM  (  /CM) 

RKt  -  RATIO  OF  MOMENTUM  *  RATIO  OF  MASS 

<*  -  A(ATOH)*  B( ELECTRON)!  G( BUFFER)!  D(ION)) 

PA*DSQRT ( 2 . 0* AMASStR VKBtTi A ) 

PB=DSQRT ( 2 , 0*EMASS*R VKBtTB ) 

PG=DSQRT ( 2 . 0*GMASS*R VKBtTG ) 

FD=DSQRT ( 2 . OtAMASStRVKBtTD ) 

RKA=1 .0 

RKB= ( PA/PB  >  t ( EMASS/AMASS ) 

RKG= <  P A/PG ) * ( GMASS/AMASS ) 

RKD= ( PA/PD )»( AMASS/AMASS ) 

RETURN 
END 

SUBROUTINE  SETTRN(NT!lS) 

IMPLICIT  REAL*8(A-HiO-Z) 

COMMON  /TRANS/  INIT r INITX » LAST r LASTX i D JE ! D JR  » 

1  GIN ! GAME ! SHIFTZ  t  NTRANS ( 50 / ! NEXCI T ( 5 ) ! FRAC i IDK 

COMMON  /ATDAT/  NSHELL ( 60 ) ! LMO  <  40 ) ! JMO  <  AO  > ! EN  <  60 ) i GAMARR ( 60 ) 


ION  11(3*50) *12(3*50) *13(3*50) *14(3*50) *15(5 *50) 

ION  R1(3.50)»R2(3.50)*R3(3>50) 

ION  NXTRAN(50) 

TINE  SETTRN  -  SET  TRANSITION  INDICES 
INPUTS 

S  -  INDEX  TO  IDENTIFY  TRANSITION 

!T  -  NUMBER  OF  THE  FINE  STRUCTURE  TRANSITION  INVOLVED 


BLES 


NTRANS(IS)  -  NUMBER  OF  TRANSITIONS  WITHOUT  EXCITATION 

NXTRAN(IS)  -  NUMBER  OF  TRANSITIONS  DUE  TO  EXCITATION 

I1(I»IS)  -  INDEX  IDENTIFYING  INITIAL  STATE 

I2(I*IS)  -  INDEX  OF  PERTURBING  FINE  STRUCTURE  LINE 

13(1* IS)  -  INDEX  OF  FINAL  STATE 

I4(I»IS)  -  INDEX  OF  PERTURBING  FINE  STRUCTURE  LINE 

15(1* IS)  -  INDEX  OF  EXCITING  STATE 

R1(I»IS)  -  DIPOLE  MATRIX  ELEMENT  SQUARED-  INITIAL 

R2(I*IS)  -  DIPOLE  MATRIX  ELEMENT  SQUARED-  FINAL 

R3(I*IS)  -  DIPOLE  MATRIX  ELEMENT  SQUARED-  6R0UND 

TRANS(18)/3/ 

XTRAN(18)/1/ 

1(1*18) >11 (2*18) *11 (3*18)/19*20»20/ 
2(1*18)*I2(2*18)*I2(3*18)/20*19*19/ 

3( 1*18) *13(2* 18) *13(3* 18)/05*05*04/ 
4(l»18)*I4(2*18)*I4(3*18)/04*04*05/ 

5( 1*18) *15(2* 18)/21*0/ 
l(l*18)*Rl(2*18)*Rl(3*18)/3t3«15E-14/ 

2(1*18) »R2(2»18)*R2(3*18)/3*1.14E-15/ 

3(1*18) >R3(2.18)*R3(3*18)/3»6.4E-16/ 

TRANS ( 24 )/3/ 

XTRAN(24)/0/ 

1( 1*24) *11 (2*24) *11(3*24 )/13*13* 12/ 

2( 1 *24 )* 12(2*24)* 12(3*24 )/12* 12* 13/ 

3( 1 >24) > 13(2*24)* 13(3*24 )/05> 04 >05/ 

4(1 *24) *14(2*24) »I4(3*24)/04 *05* 05/ 
l(l>24)»R2(2*24)»R3(2*24)/3*2.0D-14/ 

2( 1 *24 )*R2(2* 24 )*R2(3*24)/3*1« 140-15/ 

3(1 *24) »R3(2*24) *R3(3*24)/3*4.4D-14/ 

I (NT *  IS) 

(2(NT>IS) 
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NXT=NXTRAN(IS) 

DO  10  IX=lfNXT 
NEXCIT(IX)=I5(IX>IS) 

10  CONTINUE 

LAST*I3(NT»IS) 

LASTX=I4(NT»IS) 

DJE*R1(NT»IS) 

DJN=R2(NT» IS) 

GIN=R3(NT»IS) 

GAME=GAMARR( INIT) 

ZCINIT=(EN(Il(l»IS))+EN(I2(l»IS)))/2.0 
ZCLAST=(EN(I3(ltIS))4EN(I4(l»IS) ))/2*0 
SHIFTZ=EN(INIT ) -ZCINIT+ZCLAST-EN (LAST ) 

DENON=2.0*(2,0*(DFLOAT(LMO(INIT) )+l .0) ) 

FRAC=  ( DFLOAT  ( JMO  ( INI T ) ) 41 . 0 )  /DENOH 
IDK=0 

IF(LHO(INIT) .NE.LMO(INITX) )IDK=1 

RETURN 

END 

SUBROUTINE  RS  I 

IMPLICIT  REALM  (A-H»0-Z) 

COMMON  /TRANS/  INIT»INITXiLAST»LASTX*DJE»DJN»  i 

1  GIN »  GAME  >  SHIFTZ » NTRANS  ( 50 ) » NEXCIT  ( 5 ) » FRAC » IDK 

COMMON  /HERDAT/  HUT (20) rHABS(20) ? NHERH  i 

COMMON  /CSNE/  AMASS » EMASS  >  GM ASS » ZET A  r  RL J 

COMMON  /RSVAL/  Rl(20)fR2(20)  J 

COMMON  /MEDM/  RHOA  » RHOB  r  RHOG  r  RHOD •  1 

1  TArTBfTGfTD, 

2  PA»PB»PG»PD,  \ 

3  RKAtRKBiRKGiRKD  \ 

COMMON  /BCON/  RITAA»RITAB»RITAG»RITAD» 

1  EPSA  r  EPSB » EPSG » EPSD » 

2  RKSIA»RKSIB»RKSIG»RKSID» 

3  C  JA » C JB » C JG » C JD » 

4  GAMJ»G»GP»DNPJ»DJ»LJ 
DIMENSION  R1E(20)»R2E(20)»R1N(20)»R2N(20) 

SUBROITINE  RS  -  CALCULATES  LINESHAPE  PARAMETERS  R1»R2 

VARIABLES 

RlE(IH)tR2E(IH)  -  LINESHPAE  PARAMETER  INITIAL  STATE 

R1N(IH)»R2N(IH)  -  LINESHAPE  PARAMETER  FINAL  STATE 

IEXCIT  -  INDEX  FOR  EXCITED  STATES 

GAME  -  RADIATIVE  UIDTH 
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IEXCIT=0 

DO  50  IH=1»NHERM 
R1<IH>*0.0 
R2(IH)=0.0 
R1E(IH)=0.0 
R2E(IH)=0.0 
R1N(IH)=0»0 
R2N(IH)=0.0 
50  CONTINUE 
60  IEXCITSIEXCIT+1 

IF  < IEXCIT . GT . 1 ) INITX=NEXCIT ( IEXCIT-1 ) 

IF(  INITX  .EQ.  0  )  GO  TO  150 

t*S  INITIAL  STATE 

CALL  SEKINITr  INITXiDJEtGIN) 

DO  100  IH=1»NHERH 
W=HABS(IH) 

IF(W.LT.O.O)  GO  TO  100 

CALL  AS(W>RITAA'EPSA'RKSIA>RKArCJA»l'ARA'AIA) 

CALL  AS(W»RITAB»EPSB»RKSIB»RKB»CJB»2*ARB»AIB> 

CALL  AS(W»RITAG.EPSG»RKSIG»RK6»CJG»3»ARGfAIG> 

CALL  AS(W»RITAD,EPSD.RKSID»RKD»CJD»4»ARD»AID> 

Rlk  *H)=R1E(IH)+ARA+ARB+ARG+ARD 
R2E( xH)=R2E( IH)+AI A+AIB+AID+AID 

PRINT  1010iRlE<IH)rARA*ARBfARGf ARDrR2E(IH)rAIArAIB»AIGrAID 
FORMAT  (IX* lHFf 1P10E12«2) 

100  CONTINUE 
GO  TO  60 
150  CONTINUE 

»*t  FINAL  STATE 

CALL  SET(LAST»LASTX»DJN»GIN) 

DO  200  IHsl*NHERM 
W-HABS(IH) 

I F(  U  .LT.  0.0  )  GO  TO  200 

CALL  AS(W»RITAA»EPSA>RKSIArRKA'CJArlrARA>AIA) 

CALL  AS( W t  RITAB t EPSB . RKSIB r RKB r  C JB » 2 >  ARB > AI B ) 

CALL  AS<U»RITA6»EPSG»RKSIG»RKG»CJG»3»ARGiAIG> 

CALL  AS(H»RITAD»EPSD»RKSIDfRKD»CJD»4rARD»AID) 

RIN(IH) =ARA+ARB+ARG+ARD 
R2N( IH)*AIA+AIB+AID+AID 

PRINT  1020fRlN(IH)?ARArARBrARGrARDrR2N(IH);AIAr4IBrAIGrAID 
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1020  FORMAT! IX* 1H-»1P10E12*2> 

200  CONTINUE 

DO  300  IH-lfNHERN 
H«HABS(IH> 

IF(M.LT.O.O)  GO  TO  300 
R1(IH)*R1E(IH)-R1N(IH) 

R2(IH)*GAME/2*0+R2E(IH)+R2N(IH) 

300  CONTINUE 
RETURN 
END 

SUBROUTINE  SET ( IDJ* IDJX* DJT » GIN) 

IMPLICIT  REAL*8(A-HrO-Z) 

COMMON  /ATDAT/  NSHELL<60)*LM0<60>» JM0(6O)»EN<60)*GAMARR(60) 

COMMON  /MEDM/  RHOA  >  RHOB » RHOG  * RHOD  * 

1  TA*TB»T6*TD» 

2  PA#PB*PG*PD* 

3  RKA»RKB»RK6>RKD 
COMMON  /BCON/  RITAA*RITAB»RITAG*RITAD* 

1  EPSA*EPSB*EPSG*EPSD* 

2  RKSIAfRKSIB»RKSIG»RKSID» 

3  CJAtCJBrCJGfCJDr 

4  GAMJ»G>6PtDNPJ'0J»LJ 
COMMON  /CSNE/  AMASS .EMASS*GMASS.ZETA»RLJ 
COMMON  /CONST/  E2»PI*VELCfRVKB»C0NSTI 

SUBROUTINE  SET  -  SETS  CONSTANTS  THAT  DEPEND  ON  MUCH  TRANSITION  IS 
TRANSITION  IS  BEING  EVALUATED 

CALL  INPUTS 

IDJ  -  INDEX  OF  PRIMARY  STATE 
IDJX  -  INDEX  OF  PERTURB IN6  STATE 
DJJ  *  DIPOLE  MATRIX  ELEMENT  SQUARED 
GIN  -  DIPOLE  MATRIX  ELEMENT  SQUARED  -  GROUND 

VARIABLES 

RITA*  -  OFF  RESONANT  ENERGY  PARAMETER 
EPS*  -  DECAY  WIDTH  PARAMETER 
RKSI*  -  MATRIX  ELEMENT  PARAMETER 
CJ*  -  COUPLING  STRENGTH 

FACTP«<2.0*PI)t*2.5 
GAMJ=GAMARR ( IDJ) 

G*GIN 


GP=*G/30. 

DJ*DJT 

DNP J=DEL ( IDJXt IDJ) 

LJ=LMO(IDJ) 

RLJ*DFLOAT(LJ) 

C 

C  ttt  ATOM-ATOM 

C 

RITAA=AMASStDNPJtDJ/12. 

EPSA*-AMASStGAM JtD J/24 . 

RKSIA=PAtDSQRT (DJ/12« ) 

C JA= 1 08 « tE2tE2t AMASStGtGtRHOA/ ( F ACTP tD JtPA  > 

C 

C  ttt  ELECTRON  -  ATOM 

C 

R I T AB*EMASStDNP JtD J/6 ♦ 

EPSfc--EMASS*GAMJ*DJ/12 . 

RKSIB-PB*DSQRT(DJ/12. ) 

C JB«3 . tE2tE2t AMASStD JtRHOB/ ( FACTPtPA > 

C 

C  tit  BUFFER  -  ATOM 

C 

R I T AG=6MASStDNP JtD J/6 • 

EPSG*-GMASStGAMJtDJ/12. 

RKSIG*P6tDSQRT(DJ/12« ) 

CJG*108.tE2tE2tAHASStGPtGPtRH06/<FACTPtPAtDJ> 

C 

C  ttt  ION  -  ATOM 

C 

RITAD=AMASStDNPJtDJ/12. 

EPSDs-AMASStGAMJtDJ/24 . 

RKSID-PDtDSQRT (DJ/12 . ) 

C JD=3 . t E2tE2t AMASStD JtRHOD/ ( FACTPtPA ) 

RETURN 

END 

FUNCTION  DEL(N1»N2) 

IMPLICIT  REALtB(A-H»0-Z) 

COMMON  /ATDAT/  NSHELL(60)»LM0<60)» JM0(60)»EN(60)»6AMARR(60> 
COMMON  /CONST/  E2»PI»VELCrRVKBfC0NSTl 
C 

C  DEL(N1»N2)  *  BOHR  FREQUENCY  (  /CM) 

C 

C  N1?N2  -  INDICIES  OF  ATOMIC  ENERGY  LEVELS 
C 


DEL*(EN<Nl)-EN<N2))t<2*0tPI) 


o  r>  o 


END 

SUBROUTINE  AS(H»RITA,EPS»RKSI»RK»CJ»IDf AR'AI) 

IMPLICIT  REAL*8(A-H»0-Z> 

COMMON  /HER DAT/  HVT ( 20 ) *  HABS ( 20 ) i NHERM 
DIMENSION  F1<20)»F2<20) 

C 

C  SUBROUTINE  AS  -  CALCULATES  DENSITY  INDEPENDENT 
C  LINESHAPE  PARAMETER 

C 

C  CALL  INPUTS 

C 

C  W  HERMITE  ABSCISSA  VALUE 

C  RITA  -  OFF  RESONANT  ENERGY  PARAMETER 

C  EPS  -  DECAY  WIDTH  PARAMETER 

C  RKSI  -  MATRIX  ELEMENT  PARAMETER 

C  RK  RATIO  OF  MOMENTUM  *  RATIO  MASS 

C  CJ  -  COUPLING  STRENGTH 

C  ID  -  INDEX  TO  IDENTIFY  SPECIES 

C 

C  CALL  OUTPUTS 

C 

C  AR  -  REAL  PART  DENS  INDEP  LINESHAPE  PARAM 

C  AI  IMAG  PART  DENS  INDEP  LINESHAPE  PARAM 

C 

DO  100  IH*1»NHERM 
V«HABS(IH) 

VH«V+RK*U 

X=DABS(VU> 

CALL  SUMRI<O.DO>XfXfRITArEPS'RKSI'ID'0>VALI) 

CALL  SUHRI(-l.D0fl.D0fX'RITA'EPS'RKSIrID»l'VALR> 
F1<IH)«SGN<VW)*VALR 
F2(IH)*SGN<VW)*VALI 
100  CONTINUE 

CALL  HERM<FlfAR> 

CALL  HERM(F2rAI) 

ARsCJtAR/W 
AI*CJ»AI/H 
RETURN 
END 

SUBROUTINE  SUMRHYLrYUfXrRITArEPSrRKSIrlDrlRrVAL) 
IMPLICIT  REAL*8(A-H»0-Z) 

SUBROUTINE  SUMRI  -  PERFORMS  A  TRAPAZOIDAL  INTEGRATION 
-  THE  FUNCTION  INTEGRATED  IS  EVALUATED 
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IN  SUBROUTINE  NBEQ 


CALL  INPUTS 

YL  -  LOWER  LIMIT  INTEGRATION 

YU  -  UPPER  LIMIT  INTEGRATION 

X  -  ABSCISSA  VALUE 

RITA  -  OFF  RESONANT  ENERGY  PARAMETER 

EPS  -  DECAY  WIDTH  PARAMETER 

RKSI  -  MATRIX  ELEMENT  PARAMETER 

ID  -  IDENTIFY  SPECIES 

IR  -  IDENTIFY  REAL  OF  IMAGINARY 

CALL  OUTPUTS 

VAL  -  VALUE  OF  INTEGATION 

DATA  NINT/lOO/ 

DY* ( YU-YL ) /DFLOAT (NINT ) 

VAL*0*0 

Y*YL 

CALL  NBEQ<XfYiRITA»£PS»RKSIr IDrIRvBVAL) 

NB*NINT 

BL-BVAL 

DO  100  IN-l.NINT 
Y*Y+DY 

CALL  NBEQ(X»Y»RITA»EPSfRKSIrID»IR»BVAL) 

VAL*VAL+ ( BVAL+BL )  /2  •  0 
BL=BVAL 
100  CONTINUE 
VAL*VAL*DY 
RETURN 
END 

SUBROUTINE  NBEO(X»Y»RITA»EPS»RKSI»ID»IR»BVAL) 
IMPLICIT  REAL*8(A-H»Q-Z> 

SUBROUTINE  NBEQ  -  EVALUATES  B  FUNCTION 

-  THE  INTERACTION  PARAMETER 

CALL  INPUTS 

X  -  ABSCISSA  VALUE 
Y  -  CURRENT  INTEGRATION  ABSCISSA 
RITA  -  OFF  RESONANT  ENERGY  PARAMETER 
EPS  -  DECAY  WIDTH  PARAMETER 


n  !->  o  o  o  o  r>  o  o  o  o  o 


C  F.'KSI  -  MATF: I X  ELEMENT  PARAMETER 

C  ID  -  INDEXIDENTIFIES  SPECIES 

C  IR  -  INDEX  IDENTIFIES  REAL  OR  IMAGINARY  PART 

CALL  OUTPUTS 


DUAL  -  VALUE  OF  INTERACTION  PARAMETER 


BVAL-*0. 

IFdR.EG.DGO  TO  100 
FAC=Y*Y-RITA/<RKSI*RKSI ) 

1F(FAC.LT.0.)RETURN 
RT=DSQRT(  (RKSI*Y/«2-RITA) 

QF--RKSI*Y+RT 
GM=-RKSI*Y-RT 
CALL  KS<GP»©PK»ID> 

CALL  K3(0M.0MK»1D) 

DUAL* 1 . O/DSQRT <  FAC ) *  <  QP*QP*GPK+QM*QM*QHK ) 

RETURN 
IOC  CONTINUE 

CALL  KS(YfYK»ID> 

HVAL*Y» YK*I*LOG ( < Y» Y+2 . 0*RKSI TXTY+RIT A»RITA ) l*2+EFS*EPS ) 

RETURN 

END 

SUBROUTINE  KS<YtVALK»ID) 

IMPLICIT  REAL*8(A-H.0-Z) 

COMMON  /CSNE/  AMASS  t EMASS  r GMASS t ZET A » RL J 
COMMON  /TRANS/  INIT»INITX»LAST rLASTXf DJErDJN? 

1  GINjGAHE»SHIFTZ»NTRANS(50) »NEXCIT(5) »FRAC»IDK 

COMMON  /BCON/  RITAA»RITABfRITAG»RITAD» 

1  EPSA  t  EPSB » EPSG  »EF'SD» 

2  RKSIA»RKSIB»RKSIG.RKSIB» 

3  CJAfCJF.CJGfCJDr 

A  GAMJ’ G»Gr j DNF J»DJ>LJ 

SUBROUTINE  Ks  -  EVALUATES  INTERACTION  KERNAL 
CALL  INPUTS 


Y  -  ABSCISSA  VALUE 
ID  -  INDEX  TO  IDENTIFY  SPECIES 

C  CALL  OUTFUTS 

C 

C  CALL  -  VALUE  OF  INTERACTION  KERNAL 


IF*. l  .NE.O.  >FAC*l./(YBY) 
ir<0ABS<Y>.LE,FAIL>FAC=0.0 

:f<id.eq.dgo  to  too 

I?;ID.EQ.3>G0  TO  100 

***  CHARGED  SPECIES  INTERACTION 

E=GAMJ*B3GRT(DJ/6.) 

POWER  - ■ -2 . 0 * ( Rl J+3 . 0 ) f  DLOG ( 1 . 0+Y* Y > ) 
VALK=DEXP<  POWER) 

IP  •  IDK  ,E0.1>  VALK=FACH‘ALK 

RETURN 

CONTINUE 

NEUTRAL  SPECIES  INTERACTION 

IF (IB.EQ.1>H*6/DJ 
IF ( ID . EQ . 3 ) H=G/ ( 30  < ID J ) 

POWER3 ( -2 . 0* ( RL J+3 . 0 ) *DLOG  < 1 . 0+Y* Y > -6 . 0* 
*DL0G(1.Q+H*Y*Y>> 

VALK=Y**4*DEXP  <  POWER ) 

IF  < IDK . EO . 1 ) VALK=VALK*FAC 

RETURN 

END 

FUNCTION  SGN<A> 

IMPLICIT  REAL*3(A-H,0-Z) 

3GN=0.0 

IF (  A  .EO.  0  )  RETURN 
3GN*A/DABS<A) 

RETURN 

END 

SUBROUTINE  HERM(FfVAL) 

IMPLICIT  REAL*8(A-H»0-Z) 

COMMON  /HERDAT/  HUT ( 20 ) » HABS ( 20 ) f  NHERM 
DIMENSION  F(20> 

***  HERMITE  INTEGRATION 

«4*  F  -  FUNCTION  TO  BE  INTEGRATED 
m  EVALUATED  AT  ABSCISSA  VALUES  HABSIIX) 
Ml  VAL  -  INTEGRATED  VALUE 


DATA  HWT/2. 22939364550-13*4. 39934099220-10*1. 08606937070-07* 
$  7 . 80255647850-06 » 2 . 2833363601 D-04 * 3 . 24377334220-03  * 

I  2.48105208870-02*1.09017206020-01*2.86675505360-01 * 

*  4 . 6224366960D-01 * 4 . 6224366960D-01 » 

*  2 . 86675505360-01 * 1 . 0901 7206020-01 * 2 . 48105208870-02  * 

4  3.24377334220-03*2 . 28338636010-04*7 . 802556478S0-06  * 

*  1.08606937070-07*4.39934099220-10*2.22939364550-13/ 
DATA  HABS/-5 . 387480890000* -4 . 603682449500 * -3 . 944764040100* 

$  -3 . 347854567300  * -2 . 788806058400  * -2 . 254974002000  * 

t  -1.733537712100*-1.234076215300*-0. 737473728500* 

t  -0 . 245340708300  *  +0 . 245340708300  * 

4  10 . 737473728500  *  1 1 . 234076215300  *+1.738537712100* 

t  +2 . 254974002000 *  +2 . 788806058400  *  +3 . 347854567300  * 

i  +3 .9447640401 00  * + 4 . 603682449500  * +5 . 337480890000/ 

DATA  NHERM  ,'20/ 

DATA  NHER2  /!0/ 
vAL=0.0 

00  100  IX-1.NHER2 

V AL=VAL+HUT  •  I X  > * ( F< IX > +F ( NHERM+ 1-IX>  > 
too  CONTINUE 
RETURN 
END 
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